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Influence of the quantum zero-point motion of a vortex 
on the electronic spectra of s-wave superconductors 
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We compute the influence of the quantum zero-point motion of a vortex on the electronic quasi- 
particle spectra of s-wave superconductors. The vortex is assumed to be pinned by a harmonic 
potential, and its coupling to the quasiparticles is computed in the framework of BCS theory. Near 
the core of the vortex, the motion leads to a shift of spectral weight away from the chemical poten- 
tial, and thereby reduces the zero bias conductance peak; additional structure at the frequency of 
the harmonic trap is also observed. 



I. INTRODUCTION 

Recent work on 'deconfined quantum critical points'^ 
between superfluid and insulating phases has pointed out 
the central role played by the quantum fluctuations of the 
vortex excitations of the superfluid state. ^■^•"''^ However, 
this work has largely been carried out using convenient 
model Hamiltonians, with little direct connection to the 
microscopic Bardeen-Cooper-Schrieffer (BCS) theory of 
the superconducting state. It is the purpose of this pa- 
per to begin an exploration of the quantum fluctuations 
of vortices in the BCS theory. We will examine the theory 
of a single vortex in an s-wave superconductor, fluctuat- 
ing harmonically in a trapping potential. Our primary 
focus will be on the influence of this quantum motion 
on the electronic local density of states (LDOS). Mod- 
ern scanning tunneling microscopy (STM) can provide 
detailed measurements of the LDOS in the vicinity of 
a vortex, and so can, in principle, detect signatures of 
vortex fluctuations. 

Existing s-wave superconductors have a large coher- 
ence length, ^, and consequently a large effective mass 
which is expected to scale as the cross-section area of the 
vortex core ~ ^^. As a result, the vortex zero-point mo- 
tion is small, and its signatures in the LDOS are probably 
indetectable. However, should a small coherence length 
s-wave superconductor be discovered in the future, the 
results here should be useful. 

We will also develop some basic formalism for ap- 
plication to d-wave superconductors. Naturally, the 
cuprate superconductors are prime candidates for observ- 
ing vortex zero-point motion with their small values of ^, 
and also indications of proximity to a superconductor- 
insulator transition (possibly of the 'deconfined' or 
'Landau-forbidden' variety) at finite doping. However, 
the application of the formalism of the present paper to 
d-wave superconductors is far from straightforward: the 
anisotropic gap and loss of rotational invariance makes 
the numerics much more demanding, and the presence 
of gapless nodal quasiparticles leads to new physics. We 
therefore defer full discussion of the d-wave case to a sec- 
ond paper— (hereafter referred to as II). 

In a general sense, our analysis can be viewed as a 



study of the influence of the phase of the superconducting 
order parameter on the vortex LDOS. However, rather 
than integrating over the phase fluctuations explicitly, 
we are encapsulating them in a collective co-ordinate, 
the position of the vortex. 

An important limitation of the calculation presented 
here is that, in determining the coupling between the 
moving vortex and the electronic excitations, we expand 
in gradients of the pairing amplitude of the electrons. 
The applicability of such an expansion places a lower 
bound on the value of the coherence length, and so we 
are not able to properly explore the limit of extremely 
short coherence lengths. A separate computation which 
does not make this gradient expansion, while working in 
an effective low energy theory, will be presented in paper 
II. 

We will begin in Section ^ by setting up the general 
formalism which couples a moving vortex to the elec- 
tronic quasiparticles. Rather than considering a vortex 
lattice and its oscillations, we will use a simple Einstein 
model for the vortex lattice phonons, and so work with a 
single vortex moving in a harmonic potential. The elec- 
tronic self-energy corrections from the vortex motion will 
be computed in Section IIIII for an s-wave superconduc- 
tor. The rotation symmetry of this case is an important 
aide in our numerical computations, and we are able to 
obtain good numerical convergence for a significant range 
of parameters. 



II. GENERAL FORMALISM 

We will develop the general formalism for a two dimen- 
sional s- or d-wave superconductor. In principle, similar 
considerations apply to three-dimensional superconduc- 
tors, but we will not present them here. Complete nu- 
merical results for s-wave superconductors appear in the 
following section, while those for d-wave superconductors 
appear in paper II. 

To model the zero-point fluctuations of a vortex in a 
superconductor with either s- or d-wave symmetry we 
use the following Bogoliubov-de Gennes like action as a 



starting point: 
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where (with h = 1) 
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A* (r - R(r)) 



A(r-R(r)) 



(2) 

Here -0(1", t) and -0(1", t) are conjugate Grassman fields 
for electrons with mass nie, Ep — kp/2me — meVp/2 is 
the Fermi energy and A(r— R(t)) is the gap operator for 
superconducting electrons in the presence of a vortex at 
position R(t) in the plane. The vortex is allowed to move 
in imaginary time r and to account for the interaction 
between vortices in a vortex lattice we put our vortex 
in a harmonic oscillator potential V^(R) with equilibrium 
position R = 0. 

For a superconductor with s-wave symmetry the gap 
operator can be taken to be a simple scalar A(r). As- 
suming there is a vortex at the center of the plane we can 
choose a gauge such that the phase of the order parame- 
ter is equal to the polar angle of the position vector, i.e. 
(with a slight abuse of notation) A(r) = A(r)e*^. 

The d-wave case is more difficult and it is customary 
to express the gap operator in terms of a non-local or- 
der parameter with a center of mass coordinate r and 
a relative coordinate r'. Expanding in powers of r' and 
keeping only terms up to second order we can eliminate 
r' and write the gap operator with 6,^2 _y2 symmetry (see 
Appendix IX|l as 
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where {a, 5} = {ah+ha)/2. This equation is equivalent to 
a gap operator with d^y symmetry as derived by Simon 
and LeefS. For a vortex at the origin we again take A(r) = 
A(r)e'^. In Appendix fXl we will discuss issues of gauge 
invariance associated with the above expressions and the 
important comments made on this issue in Ref. |3. A 
choice of gauge is made in Eq. |(2Jl, and no additional term 
is needed to preserve gauge invariance of the Bogoliubov- 
de Gennes equations. 

Integrating out the electronic degrees of freedom in the 
action given in Eq. (^ one obtains an effective action 
for the vortex degrees of freedom. Expanding in powers 
of the velocity of the vortex we obtain for the effective 
vortex action (for simplicity at zero temperature) in the 
imaginary frequency formalism (see Appendix IB|| 
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While the term proportional to lo can be identified to be 
the Magnus force, the term proportional to uP' is just the 
vortex kinetic energy and defines the mass of the vortex, 
TO„. For a BCS-superconductor with s-wave symmetry 
va^ is of the order of rrieikp^)'^ . With ruy given, we have 
defined ujq by V{IV) = myUjQll'^ /2. For a vortex in a 
vortex lattice, the characteristic frequency is the plasma 
frequency 



ujp = \/ATr^ps/myAo 



(5) 



where ps is the superfiuid stiffness and 1/Aq is the density 
of vortices. The frequency luq is approximately given by 
u}Q w (5/2)wpi^ Finally, in a Galilean invariant superfiuid, 
LUc = 2Trps/my and in a dual picture can be identified to 
be the 'cyclotron' frequency. For a superfiuid on a lattice, 
however, it has recently been argued that lUc is reduced 
and the density of the Mott insulator has to be subtracted 
from the superfiuid stiffness when the superfiuid is close 
to a nearby Mott insulating statei^iSiii^ 

Diagrammatically, the above propagator to the effec- 
tive vortex action is just the RPA propagator which 
includes the Berry phase term. To determine the ef- 
fect of the vortex motion on the electronic spectrum we 
can use the RPA propagator and calculate the self en- 
ergy correction to the electronic eigenstates in the GW 
approximation. ^° Let us first expand the gap operator 
A (r - R(t)) in Eq. (QJ to leading order in R(r). Then 
S can be written a,s S ~ So + Sint with 
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and Tig^jQ obtained from Eq. by setting R = 0. To 
leading order the coupling of the vortex to the electronic 
degrees of freedom is described by 
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(7) 
Let us now go to a basis which diagonalizes the 
Bogoliubov-de Gennes Hamiltonian Tigdc- ^^ ^^ well- 
known that if ^f(r) = [u£(r), ^^(r)]^ (with, say, £ > 
0) is an eigenstate of Tig^Q with eigenvalue e^, then 
[— i^Kr), u|(r)]^ is an eigenstate of Hg^G '^^^^ eigenvalue 
—ei}^ Defining 
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we can write 
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such that 
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where Oz (like ax which we will use below) is just a 
usual Pauli matrix. While it is customary to restrict 
the energies tt to values greater zero^i (for which we 
use quantum numbers i > 0), we prefer to subsume 
the 'spin' index into the index i and define for £ > 
^-i(x) = {u-i{x),v-i[v)Y = {-vl(v),u\[v)Y such that 
t-i — —tg. With this notation, 






^*£(r)x£(r) 



and 5o reduces to 
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where the sum is now over all quantum numbers i in- 
cluding the 'spin' index. In the new basis, Smt can be 
written as 



S,nt = - /"dT^R(T)-Mf; 



e'XeXi' 



(13) 



where the transition matrix elements 

Mf:f> = I d^r (u*idrAve, +v^drA*uA (14) 



are like R(t) vectors in the two-dimensional plane. It is 
convenient to also define M^ = {M^ ±iM^)/2 such that 



(15) 
with dz = {dx + idy)/2. We can now calculate the self 
energy in the GW approximationiS for which we obtain 



E£(i(D) = V -:: 
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with 
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Here, a = ± and ujy = \/oJq -f ojI/A is the vortex 
('magnetoplasma') frequency in an Einstein model. As 



should be expected, the self energy satisfies Ti-i{iuo) — 

-T;t{-iLd). 

If our system is infinitely large and boundary effects 
can be neglected we can make use of the Hellmann- 
Feynman theorem 

I d\ ^\ d^nlaG *^' - i^i' - ^e) I d^r ^Idr^i, (18) 



and write M"^, as 
with 



(19) 
(20) 



In our numerical calculation presented below, however, 
we will, of course, consider a system of finite size such 
that boundary effects need to be included. We will see 
that calculating M"^, using the Hellmann-Feynman the- 
orem will turn out to be easier than calculating these 
transition matrix elements directly. 



III. VORTEX IN AN s-WAVE 
SUPERCONDUCTOR 

For a vortex in an s-wave superconductor centered 
at the origin we choose A(r) = A(r) e*^ with A(r) = 
Ao tanh(r/^) where Aq is the bulk gap and ^ ~ vp/i^Ao 
is the coherence length. 'H%,^q is rotationally invariant 
such that angular momentum is a good quantum num- 
ber. Following Caroli, de Gennes, and Matricon,^^ we 
denote angular momentum hy fi — ±1/2, ±3/2, . . . and 
write 



K{v) 

«;:(r) 



exp[-i(/z- (7z/2)6l] 



27r 






(21) 



There is only one bound state for each angular momen- 
tum /i, but there are also extended states, such that we 
also include a radial quantum number n^ It should be 
noted that while € is a collective label for jjl and n, —I 
collectively labels —^ and n. Making the above ansatz, 
the Bogoliubov-de Gennes equations reduce to 
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(22) 



Before solving these equations numerically, it is worth- 
while to first consider the bound states with energies 



e^ much smaller than the bulk gap Aq and large coher- 
ence lengths £,. As was shown by Caroli, de Gennes, and 



Matriconi*, the bound state energies are then given by 



e^ = Eifi, with Ei^-^ . (23) 

Also, the radial part of the wave function is well approx- 
imated by 



where Jm{x) are ordinary Bessel functions of the first 
kind and integer order m = fj, ^ 1/2. The constants 
C^ are independent of the ± index and are all of order 
\/kpl\. It should be noted that /^+i,_|_(r) « f^j,-{r) 
which we will use for an estimate of the matrix elements 
With the radial quantum number 



included we have 
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dr 



{[rdrA ~ A(r)] /;:,+(r)/;,;_(r) + [rdrA + A(r)] /;:,_(r)/;,;+(r)} . (25) 



r 



Using the above, we obtain yl+ « Vp/{4:Tr'^mvUJv£,'^)- 
Since the mass of a vortex in a BCS superconductor, 
my, is of the order meikpS,)'^ , it follows 



At 



A = 
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47rtii7mf(wi,/Ao)C^ 
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The fact that if we keep cj^/Aq constant, the self energy 
increases with the fifth power of 1/^ as ^ decreases is 
quite remarkable. For large ^, the self energy correction 
is very small and the motion of the vortex obviously has 
practically no influence on the spectrum. However, as S, 
decreases, the self energy correction becomes more and 
more important and we expect a dramatic change of the 
spectrum within a small range of the coherence length ^. 
Let us now consider the local density of states (LDOS) 
which (in the more general case and with (. including the 
'spin' index) is given by 



p(r,^; 






udvW 



uj — ei — 'Se{u>) + «0+ 



(27) 



For the case of an s-wave order parameter as considered 
here, the only bound state wave function which does not 
vanish at r = is m^=i/2(^)- We can therefore actually 
calculate the LDOS at the vortex center and obtain 



p{r = 0,w) 



V=l/2 



(0)P 



{{ujy +ujc + Ei)'^d{uj - £1/2) 
+ AS{uj-ei/2 - LOi) + A6{lj - ei/2 + cji)} , (28) 



u;f 



with llIi — y/icLiy + ^c + Ei)'^ + 2A which for large co- 
herence lengths ^ is close to the vortex frequency cuy. 
While there is just a peak with weight |M;i=i/2(0)P at 
the unperturbed energy ei/2 in the absence of the matrix 
element ^, as ^ decreases A increases and we also find 
two satellite peaks shifted from this position by ±uji. 

To obtain the LDOS away from the vortex center and 
to calculate the LDOS at arbitrary energy we follow Gygi 



and Schliiter^ and evaluate the eigenenergies and eigen- 
vectors of TigjQ numerically. First we replace the sys- 
tem of infinite size by a disk of finite radius Rq . We then 
expand the quasi-particle amplitudes /" _|_ into Fourier- 
Bessel series: Denoting the j'th zero of the ordinary 
Bessel function J™ by amj, we introduce the functions 



(f>mjir) 



V2 



Ro\Jm+l{oi 



rJmiamjr/Ro) , (29) 



m] ) 



which are eigenfunctions to the kinetic energy operator 
and satisfy the normalization condition 



flo 



dr r (j)mj ir)4>mj' (r) = Sjj' 



(30) 



Truncating the Fourier-Bessel series for the quasi-particle 
amplitudes /"_|_ at large Nq, we have 



/u.-M 
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With these approximations the Bogoliubov-de Gennes 
equations can be solved by solving the following matrix 
eigenvalue problem: 



1 \T/" — e'HT/" 
^T y+ j "I'm - ^M^M 



(32) 



Here T^ and A are A^o x Aq matrices with matrix ele- 
ments 
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(33) 



Ajjv = / drr ^^,-i/2,j (r) A(r) ^^,+1/2.3' {r) , (34) 




and *JJ is given by *^ = (ci • • • catq, di • • • dArJ^. Having 
calculated these eigenvectors we can express the matrix 
elements M^„.^_|.j „, which determine the A^.„„/ as 
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(35) 



where we have defined 

K^ ^ f "dr [rdrA ~ A(r)] 0^_i/2.,(r) 0^+3/2.,' (r) , 
Jo 



and 



(36) 



l}g = / dr [rdrA + A(r)] 0^+i/2,, (r) (/'^+i/2,,' (r) . 
Jo 

(37) 

Alternatively, as discussed at the end of the preceding 

section, we can also make use of the Hellmann-Feynman 

theorem which when including boundary terms reads for 

the s-wave case considered herein 



M+^, = [et ~ ei)U+f, 
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Expressing all wave functions in terms of their Fourier- 
Bessel components both integrals can be done analyti- 
cally and we obtain 
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with the matrix elements JC", and C™, given by 
^33' = '^^Sn (w + 1/2) (-!)■' ■' ^ ,.2 —■ 



^o(a™+ijv -al,j, 



(40) 
4™) = sign (m + 1/2) (-1)^-^' ""-^""t^-^" . (41) 



For simplicity and to avoid too many parameters which 
do not change the essential physics, we have neglected 
the electromagnetic vector potential in our above con- 
siderations and will not calculate the pair potential self- 
consistently. Neglecting the vector potential is safe in 
the extreme type II case where the London penetration 
length is much larger than the coherence length which 
we are considering here. Although there are significant 
deviations from the tanh-behavior of A(r) for small tem- 
peratures, corrections to the LDOS are expected to be 
small and can easily be incorporated. Once all eigenener- 
gies and eigenvectors oiTi'^^Q are determined, the matrix 
elements A^.^^, and then the LDOS can be calculated. 



In STM experiments the local tunneling conductance 
G — dl/dV can be measured as a function of gate volt- 
age V. If we are at very low temperature and have a 
tip with a constant density of states (DOS) the tunnel- 
ing conductance is essentially equal to the LDOS of the 
probe. At finite temperature, each peak in the LDOS be- 
comes broadened with the derivative of the Fermi func- 
tion f{uj) = l/(e"/^ + 1) such that 



G(r, uj = eV) = ^— j duo' p(r, oj + J) f'{uo') 
Po J 



(42) 



Here we have expressed the normalization constant in 
terms of the DOS of a free 2-dimensional electron gas 
(per spin direction), po — "^e/27r, and the corresponding 
tunneling conductance, Gq. 

In our numerical calculations we have chosen kpR — 
400. For the vortex (Einstein) frequency we use LOy — 
0.2 Aq and for simplicity we set ujc — 0. Setting the 
temperature equal to T = 0.006 Ep we obtain smooth 
curves for the tunneling conductance with well distin- 
guished bound state peaks. In Fig. ^ we show plots of 
the tunneling conductance G{r, to) at the vortex center 
and at several distances away from it for kp^ = 10 and 
rriy — mf.[kp(^Y — 100 TOg. The plots arc almost indistin- 
guishable from plots for the tunneling conductance with 
the same parameters but a static vortex (this amounts to 
taking the limit 77i„ -^ {»),- As we decrease the coher- 
ence length (and/or the vortex mass) the influence of the 
moving vortex on the electronic spectrum becomes more 
and more pronounced. In Fig. 13 we show the tunneling 
conductance for kp^ = 5 and nriy = me{kpS,)'^ — 25r7ie. 
Now one can clearly see two satellite peaks at the cen- 
ter of the vortex, shifted from the central peak by ztwi 
which is of the order of the vortex frequency LUy . In units 
of the Fermi energy the temperatures chosen in Figs. ^ 
and 13 are the same. However, in units of the bulk gap 
Aq, the temperature in the latter figure is smaller than 
that in the former by a factor of two. As a result of 
this, the hight of the zero bias peak is actually larger 
in Fig. El than in Fig. ^ even though spectral weight is 
shifted to the satellite peaks in Fig. |21 In fact, the quan- 
tum zero-point motion does lead to a reduction of the 
zero bias conductance peak as can be seen by comparing 
the results of Fig. |21 with those calculated for the same 
parameters, but with a static vortex. This case is shown 
in Fig. 13 



IV. CONCLUSIONS 

We presented a first determination of the infiuence of 
the vortex zero-point motion in a simple model based 
upon BCS theory. For s-wave superconductors, a com- 
prehensive computation was possible, with good conver- 
gence of finite-size effects, and the ability to explore a 
wide range of parameters as T ^ 0. We found two im- 
portant effects in the electronic LDOS: (z) a suppression 
in the strength of the zero-bias peak at the vortex core. 
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FIG. 1: Tunneling conductance G — dl /dV for a supercon- 
ductor with s-wave symmetry as a function of uj at the vor- 
tex center r = (upper curve) and kpr = 4, 8, . . . , 20 for 
kp^, = 10 and m„ = me{kF£,)^ = 100 me. Each curve is offset 
by 2 units for clarity. We have chosen uJv = 0.2 Ao and have 
for simplicity set uJc = 0. The temperature is T = 0.006 Ef 
and this finite temperature leads to a broadening of the con- 
ductance peaks. As can be seen in the figure, there are no 
satellite peaks. For these parameters, the tunneling conduc- 
tance of the static vortex case, without the self energy correc- 
tion, has an almost identical appearance. 



FIG. 2; Tunneling conductance G = dl/dV for a 
smaller value of ^ than that in Fig. kpS, ~ 5 and 
m„ = m^{kpi_Y ~ 25 me. The tunneling conductance is 
plotted as a function of u) at the vortex center r = Q (upper 
curve) and kpr = 4, 8, ... , 20. All other parameters are 
chosen as in FigQ The two satellite peaks shifted from the 
central peak by ±lu\ can clearly be seen in the curve for 
r = 0. 



and (ii) the appearance of satellite peaks at frequencies 
of order ±aj.u , where w^ is the vortex oscillation frequency. 

An extension of our results to d-wave superconductors 
appears in paper II. There we will introduce a model de- 
signed specifically to study the influence of the low energy 
quasiparticles, and present implications for experiments 
on the cuprate supercondutors. 
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APPENDIX A: BOGOLIUBOV-DE GENNES 

EQUATIONS FOR A d-WAVE 

SUPERCONDUCTOR 

In this appendix we derive local Bogoliubov-de Gennes 
equations for a superconductor with d-wave symmetry. 
These are applied to the vortex zero point motion, as in 
Section m while complete numerical results are deferred 
to paper II. 

As a starting point we use the following generalized 
Bogoliubov-de Gennes equations, 



Ho A 

A* -n*Q 






<^l 



ue{r) 
W£(r) 



(Al) 



Here, 



Ho = t;^ {-^dr + A{r)f - Ef (A2) 

is the effective single particle Hamiltonian, A(r) is the 
vector potential, and the gap operator A = A(r) is de- 




FIG. 3: Tunneling conductance G = dl /dV for the same 
parameters as in Fig. |5| but with a static vortex. 



fined by 

A(r)w(r) = / dV A(r - v' /2, r') v{v - r') , (A3) 

where A(r, r') is the pair potential which is a function of 
the center of mass coordinate r and the relative coordi- 
nate r'. More explicitly, A can be written as 



A(r) 



dVA(r-r72,r')e" 



-r'Sr 



(A4) 



As can easily be checked, the Bogoliubov-de Gennes 
equations IjAip are invariant under the following gauge 
transformation: 



■'x(r) 



u(r) 



u{v) — > ( 

v{r) -> e-'^('')w(r) , 

A(r,r') -> e'x(r+r'/2)+»x(r-r'/2)^(j.^j./) ^ 

A(r) -> A(r) - d,x{r^) . (A5) 

All of the above is quite general and the s-wave case is 
easily recovered by setting A(r, r') = A(r)J(r'). 

For a d-wave superconductor we express the relative 
coordinate r' in polar coordinates r' and 6' . Neglecting 
higher harmonics we write the (i^2_j,2-wave order param- 
eter as 

A(r, r') = 2A(r, r') cos(26'') . (A6) 

It is important to note that Eq. ljA6p already implies 
some choice of gauge. 



To obtain local Bogoliubov-de Gennes equations we use 
this order parameter and expand the rhs of Eq. (|A4() up 
to second order in r' which certainly is a good approx- 
imation if the dominant contribution to A(r, r') comes 
from small r' which we will assume here. Within this 
approximation, it is straightforward to show that 



{dl 



dl) 



dr^d^ 



dyAdy 



Hdl^dl) 



Akj. 



kp 



kp 



which can also be written as 



A = 



R,R,A(r)}} {9„{a„A(r)}} 



hjp 



kp 



(A7) 



(A8) 



where {a, 5} = [ab + ba)/2 denotes a symmetrized prod- 
uct and we have defined A(r) as 



/■OO 

A(r) = -4 / dr'r'^A{r,r') 
4 Jo 



(A9) 



Eq. IJA8|I is equivalent to the gap operator derived by Si- 
mon and Lee^ for a dxy-^scve superconductor. It was later 
claimed by Vafek et al.^ that this gap operator would 
not preserve the gauge invariance of the Bogoliubov-de 
Gennes equations and an additional term was introduced 
to fix this problem. The implicit assumption underlying 
this claim is that A(r) (as well as A) should transform 
under a gauge transformation as A(r) -^ e^*-^*^''-' A(r). Us- 
ing our approach, it is however easy to see that although 
the gap operator A is gauge invariant the above assump- 
tion does not hold: Working to the same order as before it 
follows from Eq. (jA5|) that under a gauge transformation 
we have A(r,r') ^ e^'^'^''^l + z9„90X?-a^V4)A(r,r'). 
Plugging this into Eq. ljA7p , it is a straightforward exer- 
cise to show that under this transformation the gap oper- 
ator A does indeed transform as A ^- e^^^^'^^ A such that 
the local Bogoliubov-de Gennes equations are gauge in- 
variant. The transformation properties of A(r) are more 
complicated but this is no problem because the gap op- 
erator is a combination of derivatives and A(r) and only 
this combination has to transform as A ^ e^*-^*-'"^A. It 
should also be noted that the derivation of the gap op- 
erator and the above argument about its gauge transfor- 
mation become exact in the limit where A(r, r') becomes 
more and more peaked at smaller and smaller r' . 



APPENDIX B: EFFECTIVE ACTION OF A 
VORTEX IN A BCS SUPERCONDUCTOR 

In this appendix we present a simple and straightfor- 
ward derivation of the effective action describing vortex 
dynamics in a clean two-dimensional BCS superconduc- 
tor. These complement the results of Ref. 0, which pre- 
sented the corresponding results using a low energy the- 
ory for a d-wave superconductor. 

Our starting point is the Bogoliubov-de Gennes-like 
action given in Eq. JQ) without the harmonic oscillator 



8 



potential V'(R), 



■01 
■01 



(Bl) 



In contrast to alternative derivations of an effective vor- 
tex action we have seen, we will not expand TiedG in 
powers of R. Instead, we diagonalize the Bogoliubov-de 
Gennes Hamiltonian TYedG at every instant in imaginary 
time r in terms of its eigenfunctions. Using the unitary 
transformation given in Eq. 111() we then obtain 



Jo f Jo f f, 



Here, Qt^i- = Qi^i'{t) is given by 



(B2) 



Q,,,. = -R(r)- / d?r [u,*(r - R(r))5rUr(r - R(t)) 

J\r\<Ra 

+ t;;(r-R(T))arWr(r-R(r))] . (B3) 

When shifting r by R(t), boundary terms which finally 
will lead to a Berry phase need to be considered carefully. 
We have therefore restricted the integration over space to 
|r| < i?o and will only at the end of our calculation take 
the limit i?o ^ oo. Doing the above mentioned shift 
r ^ r + R(t), we get Qi,i'{t) = Qfl{T) + Q''P^,{t) with 



Q("l(r) = -R. 



r|<-Ro 



d?r \uldrUi> + v}drViA , (B4) 



and 



Qtl'i.'^) ^ f {dr xR) ■ e^ (R • [u*^drUi' + v'^drVt] 

J\r\=Ro ^ 

(B5) 
To obtain an effective action for the vortex, we can now 
integrate out the fermionic degrees of freedom. It is con- 
venient to first transform the imaginary-time integral in 
Eq. IJB2|) in a Matsubara sum over fermionic frequencies 
Wn = (2n -t- 1)//?. The effective action for the vortex is 
then given by 



cVortcx^_^j.j^( J 



GQ 



(B6) 



Here, G is a diagonal matrix Green function with matrix 
elements l/(iJ)„ — Cf), Q has as its matrix elements the 
Fourier transforms of Qf ,f ' (r) , and the trace is over both 
LOn and all i. Expanding the logarithm in the effective 
vortex action up to second order we obtain 5^^'*°'^ « 

cVortcx I cVortcx ,,,ifVi 



5c'^;r = E/(^^)QS(o) 



(B7) 



and 



<? Vortex 
'-'eff,2 



11 y 



/(Q')-/(q),q(0)(^^^)|2 



££ + Cf' 



(B8) 



Here, /(e) — XjieP"^ -f 1) is the Fermi function and the 
ujra — 2TTm/ (3 are bosonic Matsubara frequencies. While 
Q^ I, does not give a contribution to Eq. HB7I) due to 

its antisymmetry, we have not included Qi\ ^,{iu}„i) into 
Eq. HB8|I because it only contributes at higher order. 
Identifying 



j(r) = : V" [uldrUi' -I- vldrViA 



2mei 



(B9) 



as the current and using well-known vector identities we 
can rewrite S, 



Vortex 
off,l 



as 



cVortcx 

'-'cff.l — 



2mei 



^ia;„ f R(iu;„i) x R(-itj,„) (J) dr ■ j 



dr(ll{-iujm)-j) 



(BIO) 



Now, using j — ps^s — Psee/^Trier, where p^ is the su- 
perfluid stiffness and v.; the superfluid velocity, we arrive 
at 

'^^'i*"'' = -^Ps^^^m iR-iiuJra) X R{-iuj^)) ■ 6^ . 

(Bll) 
This is the Magnus force whose interpretation as a ge- 
ometrical Berry phase is most transparent when trans- 
forming back to imaginary time. Defining 



1 f^ 



(Rx r") 



1 



{dR X R) • e^ 



(B12) 



as the area enclosed by the loop F we can recast Eq. (jBlip 
into 



Vortex 



oVon 
'-'eff,l 



'i2TTpsS(T) 



(B13) 



which agrees with the result by Ao and Thouless who 
emphasized the robustness of the Berry phascii^ Let us 
now consider 5^2*°'^ which we rewrite as 



cVortcx 
'-'cff,2 



1 Y^ /(£<") -/(q) 2 ITT vr M2 



2/3 



££ + ££' 



^£,1' — / d^r [uldrUgr + v}drVt 



(B14) 



(B15) 



with 



Neglecting the small iujm in the denominator of Eq. HB14|) 
and noticing that after summation over i and t' the term 
proportional to {R.{iuJm) x R{—iuJm)) ■ &z has to vanish 
we obtain 



(B16) 



ruy 1 



where the mass of the vortex, mt,, is given by 



^E 






U, 



(B17) 



Using the formahsm described in the main body of this 
paper this equation can be used to calculate the mass of a 
vortex in a superconductor with s- or d-wave symmetry. 
To summarize, we can write the vortex action as 



cVortcx 

'JnfF 



~2 






a 






, ,2 



(with ujc = 2'iTps/m^) which when taking the limit of 
zero temperature and putting the vortex in a harmonic 
oszillator potential turns into Eq. Q. While the vortex 
mass can in principle be calculated microscopically it can 
also be treated as a phenomenological constant. 
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